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We consider a simple model of cosmological collapse driven by canonical fields with exponential 
potentials. We generalise the two-field ekpyrotic collapse to consider non-orthogonal or tilted po- 
tentials and give the general condition for isocurvature field fluctuations to have a scale-invariant 
spectrum in this model. In particular we show that tilted potentials allow for a slightly red spectrum 
of perturbations as required by current observations. However a red spectrum of fluctuations implies 
that the two-field ekpyrotic phase must have a finite duration and requires a preceding phase which 
sets the initial conditions for what otherwise appears to be a fine-tuned trajectory in the phase 
space. 
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Introduction Understanding the origin of structure 
in our Universe is one of the biggest challenges in modern 
cosmology. An inflationary expansion in the very early 
universe has become the standard explanation, address- 
ing the flatness and the horizon problems as well as seed- 
ing an almost scale-invariant, nearly Gaussian distribu- 
tion of inhomogeneous perturbations about a Friedmann- 
Robertson- Walker spacetimc [1, 2]. Nonetheless, it is in- 
teresting to ask if there are alternative scenarios that 
can source primordial perturbations consistent with cur- 
rent observations. We require primordial density pertur- 
bations which are well-described by a power spectrum 
V ( (k) cx fc"^ 1 where 0.944 < n c < 0.992 [3], and the 
distribution must be sufficiently Gaussian, such that the 
amplitude of the bispectrum with respect to the square of 
the power spectrum, given by the non-linearity param- 
eter /nl, is constrained to be —10 < /nl < 74 [3] for 
local-type non-Gaussianity [4]. 

Pre-Big Bang models offer a possible alternative where 
the comoving Hubble-horizon shrinks during a collapse 
phase, generating a distribution of classical fluctuations 
on super- Hubble scales [5, 6]. One of such model is an 
ekpyrotic collapse prior to the Big Bang [7-9] where a 
canonical scalar field with a steep, negative potential en- 
ergy drives the contraction. The potential for this field 
is taken to be V(<p) = — Vbexp(— c<j>) with c 2 3> 1, which 
has a scale-invariant form, leading to a power-law col- 
lapse (a cx (—t) p where p = 2/c 2 ) and a power-law power 
spectrum of fluctuations. All collapse models face a chal- 
lenge to connect this runaway collapse to a decelerated 
expansion, but in any case this single- field model predicts 
a steep blue spectrum of adiabatic density perturbations, 
riQ ~ 3 [10], in contradiction with observations. 

An almost scale-invariant distribution of perturbations 
can be realised in the new ekpyrotic scenario [11-13] 
by considering a multi-field system [14]. Each field 
has its own steep, negative potential, e.g., V(0i,<fo) — 
— V\e" Cl ^ 1 — V2e~ C2 ^ 2 . One can perform a rotation in 
field space and define an adiabatic direction, er, and an 
isocurvature direction (x) [14-17], i.e., 

V(a, x) = -e- ca {v x e- {c ^ c ^ cx + U 2 e (c2/ci)cx ) , (1) 
with c~ 2 = ^2iC^ 2 . A power- law solution exists where 



the field \ sits m the extremum of the potential, xoi yield- 
ing a power-law solution driven by the exponential depen- 
dence, V cx e _CCT . The adiabatic fluctuations have a steep 
blue spectrum as before, but isocurvature fluctuations 
can also source the primordial density perturbation. The 
isocurvature field spectral tilt is given by n x — 1 ~ 4/c 2 
during an ekpyrotic contraction with c 2 3> 1; therefore 
the power spectrum can be nearly scale invariant. 

An essential feature of this two-field model is that the 
the power-law solution solution with p — 2 /c 2 is unstable; 
there is a tachyonic instability since the effective mass- 
squared of the x field is negative [17]. Such an insta- 
bility is necessary to achieve an almost scale-invariant 
spectrum. Quantum fluctuations on the Hubble scale 
have a power spectrum V x ~ (c 4 /4) (H/2ir) 2 which grows 
rapidly during collapse, therefore the power spectrum on 
larger scales must also experience a rapid growth, pro- 
portional to H 2 , in order to keep pace with the growing 
power on the shrinking Hubble scale. 

This raises the question of how the universe started 
sufficiently close to this unstable solution, which we will 
return to later. However the tachyonic instability does 
provide a mechanism to convert isocurvature field fluctu- 
ations into density perturbations [18]. The growth of the 
X field leads to a change from the two-field solution with 
p = 2/c 2 to a single- field solution with either p = 2/c 2 or 
p = 2/c 2 . The corresponding change in the local equa- 
tion of state, controlled by the local value of the x field, 
leads to a density perturbation, ( cx 5x [18]. Other mech- 
anisms have also been proposed which could convert the 
isocurvature field fluctuations to density perturbations 
including a kinetic conversion due to an abrupt change 
in the field trajectory after the ekpyrotic phase [11] or 
a curvaton-type conversion due to modulated reheating 
in an expanding phase following the bounce [19]. In any 
case any linear process preserves the scale dependence of 
the power spectrum and we have uq — n x . Note how- 
ever that the power spectrum is slightly blue, tiq > 1, in 
tension with current observations [3] . 

Non-linearity in the evolution of perturbations also 
provides important constraints on the model. The tachy- 
onic conversion of isocurvature field fluctuations into den- 
sity perturbations leads to local-type non-Gaussianity 
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characterised by the nonlinearity parameter [20] /nl = 
— (5/12)c 2 for 7 = 1,2. Given that we must have c 2 > c 2 
this implies /nl < (— 5/3)(l — n^) -1 , e.g., if 1 — < 0.01 
we require /nl < —100, in contradiction with obser- 
vations. Alternative conversion processes can lead to 
model-dependent results for non-Gaussianity and in par- 
ticular the kinetic conversion can lead to /nl ~ ±c which 
may be compatible with observational constraints given 
above. 

In this letter we will look at consequences of simple 
generalizations of the new ekpyrotic scenario to include 
non-orthogonal potentials and how this alters the pre- 
dicted distribution of super-Hubble perturbations and 
the problem of initial conditions. 

Tilting the potentials We will consider n canoni- 
cal scalar fields, <f> — (</>i, (j) n ), with m exponential 
potentials 



(2) 



J=i 



where c 7 = (cji, . . . , Cj n ). We recover the new ekpyrotic 
model, described above, as a special case of two orthogo- 
nal vectors c x .c 2 = 0, but in the following we will consider 
the more general case of non-orthogonal or "tilted" po- 
tentials, such that Cj.Cj ^ [14]. We restrict our discus- 
sion to Vj > so that every term in Eq. (2) is negative 
and V < 0. The case of positive potentials, V > 0, was 
discussed previously in the context of assisted inflation 
[21, 22]. We will assume that the m different vectors, 
Cj, are linearly independent. Hence our analysis is also 
restricted to m < n and we assume that the fields are 
not trapped, so that there always exists a regime with fi- 
nite energy density in which Vje~- J '- — > for any given 
potential. 

We note that we could choose to work with fields tpj cx 
Cj.cj) aligned with the potentials in Eq. (2) but then these 
fields would have a non-diagonal metric in field space, i.e., 
be non-orthogonal for Cj.c j ^ 0. 

The evolution equation for the canonical fields is given 

by 



' + ' iH 't + cj v J e ~ £j - = 



where the Friedmann equation for H = a /a is 



(3) 



(4) 



We have set 8irG = 1 and dots correspond to derivatives 
with respect to cosmic time. 

Firstly let's look at the stability of this type of system. 
To do so let's follow [17, 23] and define 



y.j = 



V6H 



je 



-cj-4> 



V3H 



(5) 
(6) 



Using Eqs (3) and (4) one finds 



dxi 
dN 

dyj 
dN 



-3x, ^1 - ^ xlj - cj.yj (7) 

VJ ^2 Xi (^ Xi ~ \ j\ Cj ^j ( 8 ) 



We can then study fixed points corresponding to scal- 
ing solutions [22-25]. They are given by dxi/dN = 
dyj/dN = 0. We find the following fixed points: 

1. Zero-potential fixed point 

These points are characterized by yj = 0,V J and 
x.x = X 2 = J2i x i — 1- This kinetic energy- 
dominated collapse with a cx (— t) 1 / 3 is unstable 
whenever there exists at least one potential with 

r 2 



- 2 >6. 



2. Single-potential fixed point 

In the case where only one potential is non-zero, 
i.e., y k 7^ while y.j — 0, V J ^ K , we then have a 
fixed point 



Vk 



K 

G 



V6 



(9) 



where (? K = J2 i c 2 Ki and we require c 2 K > 6. This 
corresponds to a power law solution of the scale fac- 
tor with a cx (— i) p , where p — 2/c|-. This collapse 
is stable with respect to the zero-potential solution. 

3. Double-potential fixed point 

In the case where two potentials are non-zero, yx ^ 
and yL ^ 0, we have a fixed point where 



Cj 

V6' 



and 



VK 



\ \°l - c k\ 



c =m\--i 



(10) 



(ii) 



with 



(c 2 K - cl_.ck) cl + (c| - c L . c K ) c K 



(12) 



\CL - CK\ 



This corresponds to a power-law solution of the 
scale factor with p — 2/c 2 where c 2 = c 2 K c 2 L — 
(c K .c L ) 2 . For this solution to exist we require 
c 2 > 6 and c K .c L < min{ cf^cl}, and in this case 
this double-potential collapse is always stable with 
respect to kinetic-dominated collapse but unstable 
with respect to single-potential collapse. 
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4. Multiple-potential fixed points 

The scaling solution for two potentials can be gen- 
eralised to the case of multiple tilted potentials. We 
again find Xi = Ci/ y/6 where we have 



c = 



(13) 




where we define the matrix Mu = Cj.Cj and 
(M -1 )/,/ is its inverse. Hence we have a power-law 
solution with p = 2/c 2 where [21] 



(14) 



A system with many exponential potentials can have 
many different fixed points, n fields with m < n poten- 
tials of the form given in Eq. (2) with independent Cj 
will have 2™ — 1 different fixed points with at least one 
non-zero potential. For instance, if we have 3 potentials 
there will be one scaling solution with 3 non-zero poten- 
tials (dk 7^ 0), three scaling solutions with 2 non-zero 
potentials, and three fixed points each with a single non- 
zero potential. In each case we can use the general result 
for the multiple-potential fixed point given in Eqs. (13) 
and (14) where the sums are to be taken over the non- 
zero potentials. This reduces to Eq.(12) for two tilted 
potentials, or c~ 2 = J2j cj 2 for multiple orthogonal po- 
tentials. 

Tilted ekpyrosis Let's consider a new ekpyrotic sce- 
nario with non-orthogonal potentials. From now on we 
will discuss the case when we have two potentials and 
two fields. Without further loss of generality we set 



ci = ci(l,0) 

£2 = c 2 (sin6>,cos#) . 



(15) 
(16) 



We recover the case of orthogonal potentials in the limit 
sin 9 — > 0. The rotation to adiabatic and isocurvature 
fields in field space is given by 



(C2 cos 9)4>i + (ci - c 2 sin 9) 4> 2 
A ' 
(ci - c 2 sin 6) (pi - (c 2 cos 9)(j> 2 



(17) 
(18) 



where A 2 = c\ — 2cic 2 sin# + c 2 . The potential in terms 
of the adiabatic field a and the isocurvature field \ IS 
V(a, x) = —e~ ca U(x) where 



c\c\ cos 2 9 



cf + c 2 ,- 2cic 2 sin# ' 



(19) 



and 



U(x) = Vi exp 

+V 2 exp 



(ci — C2 sin 9)c 
C2 cos 9 X 

(c 2 -ci sinfl)c 
ci cos 9 & 



+ 



(20) 



which reduces to (1) when 9 — 0. We note that V\,V 2 > 
so U (x) is bounded from below and has a minimum at 
X — Xo f° r cic 2 sin# < min{ cf,c 2 }. Thus there is an 
ekpyrotic, power-law solution with \ — Xo and V oc e _c<T . 
Around the minimum we can expand U(x) as 



where 



U(x)^U l + ^(x-Xo) 2 + 



/i 2 = c 2 - cxc 2 sin 6* > . 



(21) 



(22) 



Quantum fluctuations in the adiabatic field, a, lead to 
a steep blue spectrum of density perturbations, = 3 in 
the fast- roll limit [10]. A nearly scale- invariant spectrum 
of perturbations can instead originate from isocurvature 
fluctuations in the isocurvature field, x- 

Linear perturbations in the isocurvature field obey 

It is convenient to work in conformal time defined by 
dt = adr. Using the Bunch-Davies vacuum state to nor- 
malize the amplitude of quantum fluctuations at early 
times (fcr — > — 00) one finds 

aS X (r) = ^ e -^+iH-ry/ 2 H^(-kr) (24) 

where Hi 2 \—kr) is a Hankel function of the second kind 
and 



9 2c 4 -2(3 + /4)c 2 + 12^ 2 



(c 2 - 2) 2 



(25) 



The scale dependence of the power spectrum is given 
by n x — 1 = dhiV x /dln k = 3 — 2i>. To obtain an exactly 
scale-invariant spectrum we require v = 3/2 and hence 



cic 2 sin 9 



3c 2 



c? — 6 



(26) 



Note that, as expected, a scale-invariant spectrum re- 
quires that the two-potential collapse is an unstable point 
in the phase space. 

In the fast-roll limit c 2 > 1 wc have 



2 1 2(3 + m 2 )c 2 - 12^ 2 

~ i 

4 c 4 



(27) 



We will focus on the fast roll, c 2 ^> 1, and small angle 
case, <C 1. In this case 



1 4 ft ■ ClC2 
n Y — 1 ~ — I 1 



(28) 



So we find a nearly scale-invariant red spectrum if 9 < 
— 3/cic 2 . The relative tilt of the potentials in field space, 
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9, alters the effective mass of the isocurvature field en- 
abling us to obtain a slightly red spectrum in contrast to 
the case of orthogonal potentials. 

Another important constraint comes from the observed 
Gaussian distribution of primordial perturbations. Due 
to the steep exponential potentials, the isocurvature field 
fluctuations do not remain Gaussian on super-Hubblc 
scales. The resulting non-Gaussianity of the subsequent 
density perturbations depends on the conversion process. 
In the simplest case the \ field perturbations away from 
the two-potential fixed point lead to a tachyonic transi- 
tion to the single-potential solution. 

Assuming there is an instantaneous transition from 
two-potential scaling solution to single-potential solu- 
tion, when the Hubble rate reaches H = Ht then 
the integrated expansion from Hi before the transi- 



wc have 



tion to H 



-2 



after 
2 In \H T 



the transition is 
+ constant. 



N = 



given by 
Following Koyama 



c\ci cos 9 

et al [20] we can compute the non-linear growth of the 
isocurvature field and the dependence of Ht upon the 
initial field fluctuations to determine the non-linear ex- 
pansion, N, and hence the bispectrum. In the fast roll 



limit we again recover /nl 



-(5/12)c 2 < -(5/12)c 2 . 



Although c is no longer uniquely determined by the tilt 
of the power spectrum, it is nonetheless required to be 
large in the fast-roll limit. On the other hand, if we re- 
quire /nl > — 10 this demands 6 < c 2 < 24. 

Initial conditions for a red tilt Having success- 
fully obtained a red spectrum of density perturbations by 
tilting the potentials we have also created a new prob- 
lem of initial conditions for the scenario as we will now 
demonstrate. 

Although a classical solution can spend an arbitrarily 
long period of time close to the unstable double-potential 
fixed point, quantum fluctuations in the isocurvature 
field inevitably lead to a tachyonic transition within a 
finite time. The transition occurs when /x 2 (<5x 2 ) <~ 1, 
where 



V x {k)d\nk, 



(29) 



where the integral is over all wavenumbers for which we 
can treat the fluctuations as effectively classical, usually 
taken to be all super-Hubble scales. For a blue spectrum 
of perturbations the integral is dominated by the shortest 
wavelengths, i.e., the Hubble scale, and we have (5x 2 ) — 
(c 4 /4) (H/2ir) 2 . Thus, taking /i 2 ~ c 2 , the transition 
must occur when \H\ < c~ 3 , but it can begin arbitrarily 
far in the past. 

But for a red spectrum of fluctuations in the isocurva- 
ture field, the variance of the field is dominated by the 
longest wavelengths. If the new ekpyrotic phase started 
infintcly far in the past, then the variance of the field on 
the largest scales would be infinitely large leading to a 
contradiction as we require the field to be close to the 
double-potential fixed point. Thus the new ekpyrotic 
phase must have started a finite time in the past and 



c 

T 



l-n 



(30) 



where fc* is the comoving Hubble scale and fcj is the initial 
comoving Hubble scale at the beginning of the new ekpy- 
rotic phase. Assuming we have a slightly red spectrum 
with 1 — n ~ 0.01 and requiring a new ekpyrotic phase 
which lasts at least 10 e- folds, i.e., fc*/fc, > e 10 , we require 
that the transition completes when \H\ < (1 — n) 1 / 2 c -3 . 

We require a phase preceding the new ekpyrotic phase 
which sets the classical background field sufficiently close 
to the fixed point, and ensures that the isocurvature field 
has a sufficiently small variance on large scales at the 
start of the new ekpyrotic phase. There are several pos- 
sibilities, one being that the isocurvature field has a mass 
parameter that changes during the evolution, inserting 
/^ 2 (cr) in Eq.(21). This could both stabilise \ — Xo at 
early times and offers another way produce a red spec- 
trum at late times [26] . However such a potential cannot 
be realised within the context of simple exponential po- 
tentials (2), and lies outside the class of scale-invariant 
potentials with scaling solutions [24]. We expect that a 
time-dependent /i 2 (a) would lead to a running of the tilt, 
n^(k) and a scale-dependent non-Gaussianity, /nl(&)- 

Within our simple model (2) one possibility could be 
a preceding phase described by a multiple-potential scal- 
ing solution that itself would be unstable with respect 
to the two-potential solution. But the spectrum of the 
isocurvature perturbations requires a careful analysis of 
the three (or more) potential system, and could be highly 
model-dependent. An alternative preceding fixed point 
already present in our two field model is the kinetic fixed 
point with vanishing potential energy. This is an unsta- 
ble fixed point but it does describe the generic behaviour 
of the system as t — > — oo. 

The kinetic-dominated fixed point where the poten- 
tials are negligible is in fact the basis of the pre-big bang 
models proposed by Gasperini and Vcneziano [5]. It is 
well-known that a kinetic-dominated collapse leads to a 
steep blue spectrum of perturbations for any massless 
fields. Thus the isocurvature field naturally has negligi- 
ble perturbations on large scales. On the other hand a 
priori there seems no particular reason why the classi- 
cal background trajectories should approach close to the 
new ekpyrotic (double-potential) solution which is a sad- 
dle point in the phase-space [17] rather than proccding 
directly to the old ekpyrotic (single-potential) solutions 
which are the stable late-time attractors. 

Finally we note that in principle we might disregard 
the ensemble average for (x 2 ) on large scales and assume 
that simply by chance quantum fluctuations away from 
X = Xo in our local patch are unusually small. This is 
unlikely a priori but one might appeal to some anthropic 
argument that only these regions are capable of giving 
rise to observers [27]. 

Conclusions In this letter we have shown that an 
ekpyrotic collapse driven by two scalar fields with non- 
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orthogonal potentials can give a scale- invariant or slightly 
red tilted spectrum of perturbations. This is in con- 
trast to the original ekpyrotic collapse with a single field 
[7] which produces a steep blue spectrum [10], or new 
ekpyrotic collapse with two orthogonal potentials [11-13] 
which yields an almost scale- invariant, but slightly blue 
spectrum of perturbations. To obtain a slightly red spec- 
trum we fine-tune the tilt such that the angle ~ 0.01 
in the fast-roll limit. 

This two-potential solution is an unstable saddle point 
in the phase-space and a red tilted spectrum of tachyonic 
field fluctuations can therefore only exist over a finite 
range of scales. Thus the two-potential solution can only 
exist for a finite time. This is possible for a particular 
class of solutions in the phase space which must evolve 
from a kinetic-dominated initial state to approach suffi- 
cently close to the two-potential saddle point. The late- 
time attractor in the phase space is a single-potential- 
dominated collapse, i.e., the old ekpyrotic collapse [7]. 

If the tachyonic transition from two-potential to single- 
potential collapse occurs then this naturally converts the 
isocurvature field fluctuations into density perturbations. 
However this potentially leads to a large, negative non- 
Gaussianity parameter, /nl < — (5/12)c 2 , in the fast-roll 



limit, c 2 3> 1, in contradiction to the observations. 

By studying a simple two-field system we have a well- 
defined model within which we can calculate the quantum 
field perturbations about classical trajectories during a 
cosmological collapse. However it leaves unanswered the 
question of whether the required tilted potentials can be 
realised within a string theory setting, as originally en- 
visaged in the ekpyrotic scenario [7], or how the initial 
state evolves sufficiently close to an unstable saddle point 
in the phase space. This would require a preceding phase 
[26], such as is envisaged within the cyclic scenario [28]. 
In all these scenarios we still need to understand whether, 
and if so how, the universe emerges from collapse to begin 
expanding and recover the standard hot big bang. 
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